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We study an Andreev quantum dot, a quantum dot inserted in a superconducting ring, with 
several levels or conducting channels. We analyze the degeneracy of the ground state as a function 
of the phase difference and of the gate voltage and find its dependence on the Coulomb interaction 
within and between channels. We compute a (noninteger) charge of the dot region and Josephson 
current. The charge-to-phase and current-to-gate voltage sensitivities are studied. We find that, 
even in the presence of Coulomb interaction between the channels, the sensitivity increases with 
the number of channels, although it does not scale linearly as in the case with no interactions. 
The Andreev quantum dot may therefore be used as a sensitive detector of magnetic flux or as a 
Josephson transistor. 

PACS numbers: 74.78. Na, 73. 21. La, 74.45. +c 



I. INTRODUCTION 

In the middle of the last century, Josephson showed 
that a nondissipative current can flow between two 
superconductors separated by an intermediate region 
when a phase difference is applied between the two 
superconductors. This intermediate region can, in 
practice, be composed of an insulator, a normal metal, 
a constriction, etc. Recent developments of nanotech- 
nology allow us to insert quantum dots, e.g., using car- 
bon nanotubes^'^ in this region. These types of Joseph- 
son junctions are called Andreev quantum dots. ' The 
charge of the Andreev dot was shown to vary smoothly 
with the phase difference between the superconductors'^'^ 
in the absence of Coulomb interaction. 

Reference showed that the charge could be tuned 
continuously from to 2e depending on the phase dif- 
ference and gate voltage V^. In Refs. and it 
was shown that the Coulomb interaction can change the 
property of the ground state, from a usual nondegener- 
ate state with two quasiparticles to a doubly degenerate 
state with one quasiparticle. The present article is a log- 
ical continuation of the above-mentioned work and is de- 
voted to the case of several channels. These channels may 
appear due to the presence of two orbits in a single- wall 
carbon nanotubes and in multiwall nanotubes , or due 
to an inhomogeneous distribution of transverse quantized 
energy levels in metal wires. In this article, we answer 
the natural questions: How is the charge dependence on 
the phase difference affected by the presence of such chan- 
nels, what is the effect of Coulomb interaction between 
channels, and how does this dependence scale with the 
number of channels? The calculation is done in the limit 
of a large superconducting gap. 

We describe the ground state properties for differ- 
ent setups with different strengths of the regular An- 



dreev reflection (AR)^'^ and crossed Andreev reflection 
(CAR) for different types of Coulomb interaction, 
such as metal nanowires with 1 channel, single-wall car- 
bon nanotubes (SWNT) with 2 orbital channels, multi- 
wall carbon nanotubes (MWNT) with 4 channels, and 
two separate SWNTs with 2x2 channels. In this ar- 
ticle, we mainly focus on the dependencies on the tun- 
able parameters gate voltage and superconducting phase 
difference. We consider the interplay between phase- 
sensitive singlet (nondegenerate) and phase-insensitive 
doubled (doubly degenerate) states from different chan- 
nels, study the current-to-gate voltage and charge-to- 
phase sensitivities, and discuss two applications, the 
Josephson transistor and a magnetometer based on an 
Andreev quantum dot.^^ 

In Sec. II, we describe a setup based on a quantum 
dot with several levels, or alternatively an intermediate 
region composed of several single- wall nanotubes/a mul- 
tiwall nanotube. The model Hamiltonian is introduced 
in Sec. III. In Sec. IV, we represent the iterative scheme 
providing the matrix elements of the A^-channel Hamilto- 
nian. Using these results, we study the degeneracy and 




FIG. 1: (Color online) Josephson junction inserted into a 
superconducting loop which is driven by an external magnetic 
flux $ and gate voltage Vg. 
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energy of the ground state as a function of the super- 
conducting phase difference and gate voltage in Sec. V. 
Furthermore, we analyze the magnetic properties of the 
dot in Sec. VI. The charge on the dot and the Josephson 
current are calculated in Sec. VII, and their "sensitivi- 
ties" are discussed in Sec. VIII. 



II. SETUP 

As shown in Fig. 1, the setup consists of a struc- 
tured Josephson junction in the form of an Andreev 
quantum dot inserted into a superconducting loop. We 
consider a quantum dot with one or several conduct- 
ing channels and account for the Coulomb interaction 
and crossed Andreev reflection. In the Andreev dot 
without Coulomb interactions, strong normal reflection 
(electron electron, hole hole) competes with An- 
dreev reflection (electron ^ hole, hole ^ electron), and 
as a result a strong dependence of the charge and cur- 
rent on the superconducting phase difference occurs. Ev- 
ery conducting channel contributes to the charge and 
current. ^^'^^ Our Andreev dots describe a situation where 
the resonance width F is much smaller than the super- 
conducting gap A (F <C A) and the length of the normal 
part L is much shorter than the superconducting coher- 
ence length ^ (L <C 0- The superconducting phase differ- 
ence (f across the junction is controlled by the magnetic 
flux <l> passing through the ring: (p = 27r<l>/<l>o, where <l>o 
is the flux quantum <l>o = 2\e\/27Th and e = — |e| is a 
charge of one electron. 

Reahstic experimental setups might be based on a 
single- wall nanotube, ' on two or more single- wall nan- 
otubes placed in parallel, or on a multiwall nanotube. 
The case with two parallel single- wall nanotubes is illus- 
trated in Fig. 2(b): This case is interesting because the 
Coulomb interaction between nanotubes can be smaller 
than in each nanotube. If such nanotubes are separated 
by a large distance, the effects of the tubes will be addi- 
tive. The second interesting and experimentally realized 
case is the multiwall nanotube with approximately equal 
interactions in each channel and between them, as de- 
picted in Fig. 2(c). 



III. MODEL 

We consider a SWNT suspended between two super- 
conducting leads, cf. Fig. 2(a). Additional gates placed 
above the nanotube allow us to define precisely the ex- 
tent of the quantum dot and therefore allow us to mod- 
ulate its energy levels. An overall gate voltage allows 
us to apply an electric field to the entire structure and 
changes the position of the normal dot levels ei (index 
i = 1, . . . , labels the channels). The normal island(s) 
can be described as zero dimensional objects. 

The model can be described by a Hamiltonian which 
includes the dot and its internal degrees of freedom, the 



leads, and the tunnel coupling between the latter two, 
H = Hj^ -\- Hq Hr^. The first part describes the 
quantum dot with N normal levels 
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where the first term represents the interaction with the 
external gate, fii = fii^ + n^^, hia- = d\^d^^^ and and 
d^^ are electronic creation and annihilation operators for 
the i^^ level in the dot, respectively; cr =t,i, and Si is 
the i^^ energy level with respect to the Fermi level. The 
second term describes the Coulomb interaction within 
the channels (Ua) and between different channels (Uij^ 
i 7^ j). The symmetric matrix Uij is positive definite, 
and all its elements are positive. The lead Hamiltonian 
describes two BCS superconductors [with a lead index 
^ = L, R (left, right)] 
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, (2) 



with an energy dispersion in the superconducting leads 
(f/c = h^k'^ /2m — and an absolute value of the gap A in 
the bulk of the superconductors. The electron hopping 
term between dots and leads is given by 



^T = ^(^l,fc7^^i + H.c.), d: 
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where 71, r = ti^,KO'ze^^^'''^^^ and t£ are tunneling am- 
plitudes between the superconductors and the dot. The 
superconductor has a phase (p/2 on the left and —(p/2 
on the right. The calculation of observables for a ther- 
mal equilibrium system starts from the evaluation of the 
partition function Z = Tr{e~^^}, where f3 is the inverse 
temperature. 

We work in the limit of a large superconducting gap 
l^^l, Uij^ Tij <C A (the so-called A ^ oo limit), where F^j 
are the resonance half- widths originating from tunneling 
processes [Eq. (3)]. In this limit, one can integrate over 
the lead degrees of freedom and benefit from the absence 
of retardation effects due to the latter. The Hamiltonian 
H can be rewritten in a simpler form 
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describing the dot alone with the superconductors defin- 
ing the boundary conditions. ^^"^'^ We expHcitly indi- 
cate the number of channels N in the Hamiltonian 
[Eq. (4)] and build an iterative scheme for its matrix el- 
ements on N. Upon integration, Eqs. (2) and (3) 
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FIG. 2: (Color online) Structure of the Andreev quantum 
dot with N channels, (a) A dot with N normal levels with 
corresponding energies Si, . . . ,Sn is inserted between two su- 
perconductors with phase difference ip across the tunnel junc- 
tions tL and tR. (b) Quantum dot system based on two paral- 
lel single- wall nanotubes. The barriers tL and tR are formed 
by two external gates. The Coulomb interaction u between 
the channels may be smaller than the interaction U inside a 
given channel, (c) Quantum dot based on a multiwall nan- 
otube. The Coulomb energy between orbital channels is ap- 
proximately equal to the ones inside each channel. 



generate the last term in Eq. (4), which implies that this 
Hamiltonian does not conserve the number of electrons. 
The coefficients Tij are directly derived from the tunnel- 
ing amplitudes ^l^r and describe two-particle tunneling 
processes involving both Andreev and normal scattering 
events at the left and right boundaries. For symmetric 
boundaries with equal transparencies I^Lp = \t^\'^ = 
we find Tij = Tij cos{(p/2). The symmetric matrix Tij de- 
scribes Andreev reflection, the annihilation of two quasi- 
particles from the i^^ and j^^ channel di^dj^ (with simul- 
taneous creation of a Cooper pair in the superconductor) 
and the opposite process. The diagonal terms Ta cor- 
respond to the regular Andreev reflection inside the i*^ 
channel, the off-diagonal terms F^j, {i ^ j) describes the 
crossed Andreev reflection (CAR) between the i^^ and j^^ 
channel. The element Tij can be treated as the resonance 
half-width of the normal tunneling process between the 
i^^ and j^^ channel. Note that for a diagonal matrix Tij 
and Uij (with Tij = 0, Uij = 0, i 7^ j) the channels are 
additive and the answers for all physical quantities can 
be obtained by summation over all channels. 



IV. ITERATIVE SCHEME 

The Hamiltonian Eq. (4) for a single channel (the first 
channel in the iterative scheme presented below) has the 
following form 



Hi = eihi + Uihl + Tildi^di^ + H.c], 



(5) 



where the nf^ contributions from the sum nf = (ni^ + 
fiii)'^ has been absorbed in a shift of ei. We denote 
the Hamiltonian of the i^^ channel as Hi and the total 
Hamiltonian of N channels as Hn- The corresponding 
matrix of the Hi Hamiltonian has dimensions 4x4 and 
can be calculated in a basis of four states = {|0)i. 



|2)i}, the state with no electrons |0), with one 
electron with spin up W) = (i||0) or spin down |^) = (i||0), 
and the two electron state |2) = (i|(i||0) 

Hi = Sim + Uinj + adiagi{fi, 0, 0, f 1} 

fi 

£1 + /7i 

ei^Ui 

fi 2£i+4/7i 
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where 



ni = ni^ + mi = Qi/e = diagi{0, 1, 1, 2} (7) 



are matrix elements of the dimensionless charge oper- 
ator and Qi = efii. In what follows, the expression 
adiagjxi, X2, . . . , x^v} stands for the N x N matrix Aij, 



where A 



and Aij = for j N - 



1. The 



eigenvalues of Eq. (6) can be easily calculated. They 
consist of two nondegenerate levels with energies 



'0/2 — ^1 + 
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and a doubly degenerate level with energy 



(8) 
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which can be split into two separate levels by an exter- 
nal magnetic field through the Zeeman effect. Usually, 
the ground state of the Hamiltonian Eq. (6) is given by 
the singlet Eq. But in the case of nonzero Coulomb in- 
teraction Ui > 0, there exists a range of parameters (p 
and Ei for which the ground state is doubly degenerate 
with energy Ei. The doublet region is defined by the 
inequality 



(£l+2[/i)2+f2<[/2, 



(10) 



where e = ei -\- 2Ui. For an asymmetric dot (I^lI 7^ |^r|), 
in the presence of electron-phonon interaction or other 
perturbations, the "critical" Ui is different from zero, 
^i,c > 0.^^ In this article, we ignore such nonideali- 
ties and concentrate on the simplest case. The same 
results are valid for any other i > 1 channel; they can 
be obtained by the substitution Si e^, Fi F^, and 
Ui Ui. 

The structure of the ground state in the one-channel 
case is presented in Fig. 3 as a function of phase (p and 
the "central" dot-level e = Si -\- 2Ui. The blue figures in 
the center represent the doublet state 9, their shapes are 
defined by the inequality (10), the white space around 
corresponds to a usual singlet state with energy (8). For 
two completely separated single-channel dots with the 
same levels Si = £2 connected in parallel, the dependence 
is the same, however, now the central figure corresponds 
to a four-fold degenerate situation. 

We now go back to the A^-channel case. We calcu- 
late the matrix elements of the Hamiltonian Eq. (4) in 
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where we keep in mind that f 12 = f 21- The matrices 



FIG. 3: (Color online) System degeneracy Vi in the ground 
state in {'^^s) space for the channel N = 1. 



the basis (g) \v)2 ^ ...^ \v) ^ = {|0,. . . ,0), |0, . . . ,t), 
|0,. ..,;),..., |2,. ..,;), |2,...,2)}, where = {|0)„ 
W)i^ II) \2)i} are states of the Hamiltonian for the sin- 
gle channel case [Eq. (6)] 
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Here, the first sum takes into account the interactions 
inside each channel (given by the direct tensor product 
of HiS and unity operators Ij in the other subspaces). 
The second sum is the interaction between channels; each 
term represents the tensor product of the charge opera- 
tor in the i^^ channel and in the j*^ channel [as given 
by Eq. (7)]. Here, 1^ is the unit matrix in the Hilbert 
subspace of the i^^ channel, and Qi is the matrix charac- 
terizing the i^^ channel charge. The third sum describes 
CAR; each term corresponds to the CAR in the i^^ and 
j*^ channels. The Coulomb interaction inside each chan- 
nel and the regular Andreev reflection are accounted for 
in the first term. 

Let us construct a recursive procedure in N for the 
matrix elements corresponding to the A^-channel Hamil- 
tonian. As the initial step of the recursion N = 1 we 
take the matrix (6). At the second step = 2 we build 
a matrix in \ v)2 in the following way: 



% = diag2{i^i,i^i +£2 + t/2 +U2i,i^i + S2 
i^i +2^2+ 4/72 +U22} 
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and the corresponding Hermitian conjugates d\^^ d\^ are 
the matrix elements of the operators d-^, d-^ and d\^^ d\^ 
in the i^^ channel basis. This is a 16 x 16 matrix (the 
Hilbert space of the 1^* channel with dimension 4 should 
be multiplied by the Hilbert space of the 2"^ one with 
the same dimension). Each "element" is a 4 x 4 block, 
"diag" and "adiag" to be in the 2"^ channel subspace. Hi 
originates from Eq. (6), with scalar elements to be mul- 
tiplied by ll. The matrices U21 = diag;^!^^ ^^12} 
and U22 = diag]^{0, 8t/i2, 8[/i2, 0} are based on the term 
U12 which is responsible for the interaction between chan- 
nels. For the noninter acting case U12 = 0, the structure 
of the ground state is presented in Fig. 4(a) [£2 — £1 = 0] 
and 4(d) [£2 — £1 = F] in coordinates ((/:', s), where 
£ = (£i+£:2)/2 + /7i + /72 is the average and renormalized 
normal level in the dot. The interacting case is shown in 
Figs. 4(g) and 4(j), where £ = (£i+£2)/2 + /7i + /72+2/7i2; 
see Sec. V for more details. 

In the same way we can write a recursive formula 
for Hat. The easiest way to write it down is 
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'ijUifij + tij [di^dj^ + 4 Ail • (13) 



Here each matrix with index i should be written as a 
matrix of size 4* x 4* (i.e., in the subspace of the i^^ 
channel). Herein each element of the matrix associated 
with the i^^ channel is multiplied by a unitary matrix of 
size 4*~^ X 4*~^ and starts on column (c — 1) 2*~^ + 1 and 
on row (r — 1) 2*~^ + 1 of the total matrix, where c = 
1, . . . , 4 is the column and r = 1, . . . , 4 is the row in the 
original 4x4 matrix of the i^^ channel. Schematically the 
recursive procedure for building the A/'-channel matrix 
can be represented as (for N = 3) 



The smallest "rectangle" corresponds to the Hamiltonian 
Hi. In a second step, we insert such blocks in the Hamil- 
tonian H2 and continue the procedure till N. At each 
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FIG. 4: (Color online) System degeneracy V2 in the ground state in the {(p,s) space for N = 2 channels. The Coulomb matrix 
Uij is defined by Eq. (16), CAR is defined by Eq. (17). The Coulomb interaction inside each channel is equal and equal to 
the resonance width, U = F. The Coulomb interaction between channels varies from left [i^ = in (a)-(c)] to right [u = U 
in (j)-(l)]; The CAR parameter varies from top [7 = in (a), (d), (g), and (j)] to bottom [7 = F in (c), (f), (i), and (1)]. 
The difference in dot levels ^£12 = £2 — £1 = F in (d)-(f) and ^£12 = in all other plots, (a) In the noninteracting case the 
degenerate region with V2 — 2 (green) inside the nondegenerate region V2 — ^ (white) is defined by the inequality (10) and 
coincides with the same region of each channel, (b) and (c) CAR smears the border and V2 goes from to 2 continuously, 
(c) Two doubly-degenerate regions (blue) shifted by ^£12; the intersection represents a fourth- fold degenerate region (green), 
(e) and (f) The border is (nonuniformly) smeared due to CAR, starting from the point where three regions with 7^2 = 0, 1, 
and 2 join together and ending with a smeared four- fold region, (g)-(i) The average Coulomb interaction between the channels 
u — lJ 12 separates apart the region with 7^2 = 1 in (a) by the distance 4i^, but in a different manner from (b); CAR smears the 
borders of the V2 = 2 region and diminishes the 7^2 = 1 regions, (j)-(l) The maximal Coulomb interaction between channels 
u = U totally splits the V2 = 2 region into two regions with 7^2 = 1 (the distance is still 4n); CAR slowly increases V2 from 
to 1 and slightly reduces the 7^2 = 1 regions. 



i^^ step, we multiply the current space by the Hilbert 
subspace of the i^^ Andreev level. 

Note that the same procedure can be realized by using 
the Bogoliubov-de Gennes equations and the scattering 
matrix approach in the tunneling regime (with a trans- 
parency of each junction much less than unity) as was 
done in Ref. .u for a single channel. However, this ap- 
proach is poorly scalable with and the calculation is 
much more cumbersome. See also Ref. 24 for the cor- 
respondence between the tunneling Hamiltonian method 
and the scattering matrix approach. 

In the present notation, the first, the second, and the 
fourth sum in Eq. (13) are diagonal. In the absence of 
superconductivity they define energy levels of the sys- 
tem, and the Coulomb interaction leads to level "repul- 
sion." The third term of Eq. (13) includes the depen- 
dence on the superconducting phase difference and mixes 
the states due to the presence of Gooper pairs in the su- 



perconductors. At = TT the Hamiltonian Eq. (13) is 
diagonal and its eigenvalues can be found analytically. 
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FIG. 5: (Color online) The Crossed Andreev reflection be- 
tween channels 1 and 2 in the fork (or Y) geometry is sup- 
pressed if the length between the NS interfaces is much larger 
then the coherence length ^. 
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FIG. 6: (Color online) System degeneracy V3 of the ground state for N = 3 channels; the same as in Figs. 4(a)-4(c), 4(g)-4(i) 
and4(j)-4(l). 



V. ENERGY LEVELS AND DEGENERACIES 

In the single channel case, the ground state is formed 
by nondegenerate [Eq. (8)] or doubly degenerate [Eq. (9)] 
levels. However, for a large number of channels one can 
expect to observe 2^^ degeneracy, where Vn is the num- 
ber of Andreev levels in the doublet state. 

To characterize the state of each channel let us first 
determine the operator 

p, = diag,{0, 1,1,0} (14) 

which is defined in the Hilbert space of the i^^ channel. 
For the case with no interaction, this operator defines the 
Pi value equal to if the i^^ channel is in the singlet state 
and 1 if the i^^ channel is in the doublet state. In the 
case including interaction, pi may vary from to 1. The 
Vn value of the operator 

N 

Vn = ^li(8)...(8)Pi0...0lA, (15) 

i=l 

characterizes the state of the entire system, ^ Vn ^ 
A^. In the noninter acting case Vn gives the number of 
channels in doublet states and can be considered as the 
degeneracy of the ground state. 

In addition, the mechanism for intrastate transitions 
with different Vn should be discussed. As mentioned 
above these transitions qualitatively change the prop- 
erties of the system, the charging effects, the trans- 



port properties, and the magnetic response to the ex- 
ternal field. Such transitions occur, for instance, due 
to electron-phonon interactions involving a continuous 
spectrum.' Their rate is suppressed by the exponential 
factor e~^/^^®, where /cg and G are the Boltzmann 
constant and temperature, respectively.^'^ Thus, transi- 
tions where the system's spin changes by 1/2 are more 
rare than transitions where the spin remains unchanged. 
Therefore, the above picture for conversion between sin- 
glet and doublet in one channel is valid for an adiabati- 
cally slow change of the parameters (p and e. 



Let us consider cases defined by different matrices Uij 
and Tij. The diagonal elements of the Coulomb ma- 
trix Uij define the interaction in each channel. The 
off-diagonal elements describing the interactions between 
channels can be of the order of the diagonal elements for 
the multichannel wire and about zero for a few separated 
wires. We can vary the "ratio" between off-diagonal to 
diagonal elements from to 1. In the same way, we can 
vary the off-diagonal elements of Tij , which are of the or- 
der of the diagonal elements and about zero for separated 
channels, e.g., see the superconductor in fork geometry 
in Fig. 5. 
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FIG. 7: (Color online) System degeneracy V4 of the ground state for N = 4 channels, (a)-(c) V4 goes from to 4 in a 
step-like manner for 7 = and smoothly for 7 > [compare with Figs. 4(a)-4(c) and 6(a)-6(c)]; the Coulomb matrix Uij is 
defined by Eq. (16) and the AR Fij is defined by Eq. (17). (d)-(f) The average Coulomb interaction between channels u = U/2 
transforms the 7^4 = 4 region to 7^4 = 1, 2, 3, and 4 regions [compare with Figs. 4(g)-4(i) and 6(d)-6(f)]. (g)-(i) The maximal 
Coulomb interaction between channels u = U splits the 7^4 = 4 region to 4 regions with 7^4 = 1 at a distance 4u [compare with 
Figs. 4(j)-4(l) and 6(g)-6(i)]. (j) Two parallel SWNT: Uij is defined by Eq. (19), Fij by Eq. (17) with 7 = 0. The structure is 



similar to the double structure in Fig. 

7 = r. 



4(g). (k) and (1) Two SWNT: Uij defined by Eq. (19), Fij by Eq. (20) for 7 = and 



A. Toy model 

We parametrize the Coulomb interaction matrix with 
two parameters 

Here U is the interaction within each channel and u is the 
interaction between any pair of channels, with u ^ U. 
The Coulomb interaction Eq. (16) describes N parallel 
wires, each with one channel (no orbital degeneracy). In 
this model the coupling to the superconductor is given 
by the matrix 

The diagonal elements F correspond to the regular An- 
dreev reflection, and off-diagonal elements 7 describe 
crossed Andreev reflection (CAR). The situation 7 < F 
is possible, e.g., in the fork geometry shown in Fig. 5. 



The results of the numerical calculation for £1 = £2 = 
. . . = En and U = F are presented in Figs. 4(a)-4(c) 
and 4(g)-4(l) [N = 2] as well as Figs. 6 [TV = 3] and 
Figs. 7(a)-7(i). Here £ = £1 + 2/7 + 2{N - l)u. 

In the absence of Coulomb interaction between the 
channels, = 0, and zero CAR 7 = 0, the channel 
contributions are additive and the system is either in 
the nondegenerate state Vn = or in the state with 
Vn = N and degeneracy 2^^, see Figs. 4(a), 6(a), and 
7(a). A nonzero 7 leads to smearing of the borders be- 
tween Vn = and Vn = N regions, see Figs. 4(b), 6(b), 
and 7(b) for 7 = O.IF and Figs. 4(c), 6(c), and 7(c) for 

7 = r. 

For < u < U and 7 = 0, the areas with integer 
Vn = 0,1, N are shown in Figs. 4(g), 6(d), and 7(d). 
The "centers" of the regions with Vn > are separated 
by distances defined by the off-diagonal elements of Uij, 
e.g., for Eq. (16) these distances are 4ii. Nonzero values 
of 7 result in continuous Vn in [0 . . . A^], see Figs. 4(h)- 
4(i), 6(e)-6(f), and 7(e)-7(f). 

For u = U and 7 = there are N regions separated 
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by 4u with Vr, = 1 [Figs. 4(j), 6(g), and 7(g)]. CAR 
7 > increases Vn in regions which were zero initiahy 
at 7 = 0; Vn can be increased as well as decreased from 
1 [Figs. 4(k)-4(l), 6(h)-6(i), and 7(h)-7(i)]. 



B. Nanotubes 

In this section, we describe Andreev quantum dots 
based on a single- wall nanotube, two parallel single- wall 
nanotubes [Fig. 2(b)], or a multiwall nanotube/molecule 
[Fig. 2(c)]. 

The single-wall nanotube specifies a two-channel 
Coulomb matrix and an Andreev refiection matrix 



'u 


u' 




"r 


r" 


u 


u 


r — 


r 


r 



(18) 



The existence of two orbital states leads to the appear- 
ance of two doubly degenerate regions, see Fig. 4(1). 

Two parallel single- wall nanotubes can be described by 
the Coulomb matrix 



and an AR matrix 



r 



U U u u 

U U u u 

u u U U 

u u U U 



r r 7 7 

r r 7 7 

7 7 r r 

7 7 r r 



(19) 



(20) 



The 2x2 block structure appears due to the twofold 
orbital degeneracy in the single-wall nanotube, where u 
and 7 describe the interaction between nanotubes. The 
structure of the ground state for 7 = is presented in 
Fig. 7(j): both regions with Vn = 1 (two- fold degener- 
acy) and Vn = 2 (four- fold degeneracy) are present. Two 
"copies" of 1-2-1 regions appear due to orbital degener- 
acy in each nanotube, which can be compared to Fig. 4(g) 
for the case with no orbital degeneracy. In Figs. 7(k) and 
7(1) the Vn behavior at 7 = O.IF and 7 = F is shown. 

The multiwall nanotube has four or more orbital states. 
The most simple case. A/" = 4, can be described by the 
interaction matrix Eq. (19) with u = U [or the same 
Eq. (16)]. This situation is illustrated in Figs. 7(g)-7(i). 

In the case of a diagonal Tij {Tij = 0, i 7^ j), the 
ground state is created from the singlet |0), |2) or dou- 
blet It), \D states (but not a mixture of singlet and 
doublet states of the same channel) and Vn takes on in- 
teger values only. In the case of arbitrary off-diagonal 
elements Tij {Tij ^ 0, i 7^ j), the ground state can be 
a mixture of all possible states including simultaneously 
|0), It), It), and 1 2) states from the same channel. 

The value of Vn has a large impact on the physical 
properties of the dot, such as magnetic properties, charge 
on the dot, and the Josephson current. 



VI. MAGNETIC PROPERTIES 

In this section, we describe the magnetic properties 
of the ground state in the presence of a weak external 
magnetic field. We take into account only the Zeeman 
splitting and neglect all other lower-order effects such as 
Rashba spin-orbit coupling (see Ref. i for more details). 
This situation can be described by an additional term 
to the Hamiltonian (4) with a magnetic field B in the 
transverse direction 



i=l cr 



(21) 



where g ^ 2 \^ the Lande factor and ji^ = \e\h/2m the 
Bohr magneton. The spin quantum number is s^/^ = 
in the corresponding i^^ channel. In the presence 
of a magnetic field, the doublet state of each channel 
splits into two levels with energy difference gii^B. 

In the absence of CAR, Vn takes on only the inte- 
ger values from to A/'. Then, the lowest level of the 
system with spin — P^/2 (where Vn is the number of 
channels in this doublet state) is nondegenerate, the first 
excited with spin —{Vn — 2)/2 is T^^r-degenerate, etc. 
The whole picture is presented in Table I. In this table 

= VnI/{Vn - i)l il and naturally EfiJ) ^^n = 2^^- 
The magnetic properties of the system correspond to a 



TABLE I: The spins and degeneracies of the states obtained 
from a 2^^ -degenerate ground state with Zeeman splitting. 



State 




Spin 


Degeneracy 


n... 
n... 


n 
n 


Vn/2 
Vn/2 - 1 


^Vn 


n... 


u 


-VN/2 + i 


^Vn 


n... 

u... 


u 
u 


-Vn/2 + 1 

-Vn/2 


^Vn ^ 



system of Vn independent spins 1/2. The infiuence of the 
Hamiltonian (21) leads to a change in size of the degener- 
ate regions, which is equivalent to changing the Coulomb 
interaction U. The dependence of Vn on (p and e has 
been discussed in Sec. V. 



VII. CHARGE AND JOSEPHSON CURRENT 

The operator for the total dot charge is given by 

N 



(22) 



=1 



-2 



2 channels w = 0, 7 ^^^1 


u = t//2, 7 = 
^^^^^^ 
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FIG. 8: (Color online) Charge of the ground state 62,3,4 in {(p,e) space. Plots (a)-(f) [N = 2], (g)-(l) [N = 3] are for Uij 
and Tij as defined by Eqs. (16) and (17). Plot (m) corresponds to the two nanotubes [Uij defined by Eq. (19)] and zero CAR 
between all channels [Eq. (17)]. Plots (n) and (o) describe two nanotubes [Uij and Fij are given by Eqs. (19) and (20)] with 
suppressed and maximal CAR between nanotubes. In all plots the charges goes from (cold red) to 2A^e (warm red); the 
"centralized" charge Ne is shown in white color. Note that for all different N there are different colorbars. 



where Qi = ehi = e {dl^d-^ + ^li^ii) chan- 
nel charge operator. The latter has matrix elements 
Qi = e diag^jO, 1, 1, 2}. The dot charge (22) has matrix 
elements 



N 

Q^ = e^n,, (23) 



where the matrices rii are inserted using the same rules 
as in formula (13). 

The charge of the ground state is presented in Fig. 8 
as a function of cp and e. In each channel, the charge 
can vary from 2e to when increasing the Coulomb 
interaction u and crossed Andreev reflections 7 change 
the total charge from the straightforward sum of channel 
charges. 
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FIG. 9: (Color online) The Josephson current of the ground state ^2,3,4 in (cp, s) space. All parameters are the same as in 
Fig. 8. A zero current is depicted by white color, negative by blue, and positive by red. For all = 2, 3,4 there are different 
colorbars. 



In the single channel case A' = 1, the charge Qi can 
vary from to 2e; for the diagonal U and F matrices 
[Eqs. (16) and (17)], the total charge Q2,3 is reduced to 
the sum of charges from each channel as shown for A" = 2 
in Fig. 8(a) and for AT = 3 in Fig. 8(g). The white flat 
plateau in the middle part of the plot corresponds to the 
region with Vn = N and = eA"; in this region, the 
charge does not depend on (p or e. ' The nonzero u "di- 
vides" this region into pieces: partially [e.g. Figs. 8(c) 
and 8(i)] or fully [e.g. Figs. 8(e) and 8(k)]. In the regions 



with Vn < only charges from A" — Vn channels con- 
tribute to the phase-dependent part of the total charge; 
the remaining channels {Vn ) give a constant contribution 
e. As previously, CAR smears the borders of the plateau 
as shown in Figs. 8(b)-8(f) and 8(h)-8(l). In Figs. 8(m)- 
8(o) the total charge Q4 of the two SWNTs connected 
in parallel is presented; for comparison see Figs. 7(j)-7(l) 
with V4. The charge in each channel is changed from to 
2e, so the total charge goes from to 2eA'. The scaling of 
the phase-dependent part of the charge will be analyzed 
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in Sec. VIII. 

The total Josephson current operator is given by the 
expression 

N 

Xa. = ^ li (8) . . . (8) /i . . . (8) Iat, (24) 

2 = 1 

where 

/, = -| r, sin I {dj^^ + dl^dl) (25) 

is the i^^ channel current operator. The latter has matrix 
elements 

Ii = -^T, sin I adiag,{l, 0, 0, 1}. (26) 

The total current operator can be generated in the same 
way as in Eqs. (13) and (23). Note that the charge 
and current can be obtained as a derivate with respect 
to and respectively. The Josephson current X^ 
corresponding to the operator (24) is presented in Fig. 9 
as a function of the superconducting phase difference (p 
and e. 

In the noninteracting case u = and 7 = 0, the cur- 
rent is defined by energy the levels [Eqs. (8) and (9)] 
and scales proportional to the number of channels see 
Fig. 9(a) for AT = 2 and in Fig. 9(g) for N = 3. The white 
region corresponds to the zero-current region; along the 
horizontal lines the current demonstrates the usual sine 
behavior (if the channel is in the singlet state) or zero cur- 
rent behavior (if the channel is in the doublet state). ^'-^^ 
In the maximally degenerate regions Vn = the current 
can fiow in the opposite direction (7r-junction) if the con- 
tinuous spectrum above A is taken into account. ^'^^"^^ 
The channels can be in the doublet state simultaneously 
[Figs. 9(a) and 9(g)] or asynchronously [Figs. 9(c), 9(e), 
9(i), and 9(k)]. The small CAR 7 smears the borders of 
the degenerate regions [Figs. 9(b), 9(d), and 9(f)]; large 
CAR decreases or even eliminates the degenerate regions 
[Figs. 9(h), 9(j), and 9(1)]. 

The total current through two parallel SWNTs is 
shown in Figs. 9(m)-9(o); for comparison see Figs. 7(j)- 
7(1). 

The Coulomb interaction as usually suppresses the cur- 
rent; CAR adds additional N{N — 1) channels for the 
current and therefore increases the current. Note that in 
the limits considered here, the maximal Vn = N leads to 
zero charge and current response with respect to (p; with 
decreasing Vn the response increases. 



VIII. SENSITIVITY 

Let us now analyze the sensitivity of the charge Q 
with the superconducting phase difference cp. The most 
relevant physical quantity is the differential sensitivity 
S = {Ae/h)d^pQ. It characterizes the charge response to 



infinitely small deviations of the superconducting phase 
and describes the operating component in new types of 
magnetometers based on the charge of the Andreev quan- 
tum dot proposed in Refs. lO and 22. The above defined 
sensitivity coincides with the current-to-gate voltage sen- 
sitivity S = edX/de^ which characterizes the Josephson 
transistor. This correspondence can easily be proven 
if we remember that the charge is the energy derivative 
with respect to the gate voltage, and the current is the 
energy derivative with respect to the phase. This implies 
that the results we find for the differential sensitivity of 
the magnetometer are also true for the sensitivity of the 
Josephson transistor. We will denote Sn as the sensitiv- 
ity of the system with TV channels. 

The sensitivity as a function of the phase difference if 
and the gate voltage e is displayed in Fig. 10. Given the 
case of multiple channels and the absence of Coulomb 
interaction between them, it is natural to expect the sen- 
sitivity to be composed from the sensitivity of each chan- 
nel separately as shown in Fig. 10(a) and in Fig. 10(g). 
In this section we examine the infiuence of the Coulomb 
interaction between channels and CAR on the sensitivity. 

Due to the slow singlet-doublet transitions, we analyze 
only the sensitivity which comes from the channels in the 
singlet regime and omit the sensitivity given by the step 
in the charge during the singlet-doublet transition: For 
low-frequency measurements this kind of sensitivity may 
be quite important. 

In the absence of the Coulomb interaction u = U = {) 
and the same normal dot levels Si = £2 = - - - = Sn-, the 
differential sensitivity has a maximum near point Lp = it 
and £ = £1 0.'" 

In case of a single channel N = 1 with Coulomb in- 
teraction U = r, the charge at = and ip = is 
insensitive to the phase. The maximal sensitivity ap- 
pears at the border of the doublet region. See the shape 
Fig. 10(a) for AT = 2 and Fig. 10(g) for N = 3. 

If we consider the two-channel case N = 2 with differ- 
ent dot levels Si ^ 82 without the Coulomb interaction 
u = U = 0, then the sensitivity is defined by the dis- 
tance fei2 = £2—Si. If this distance equals to zero, then 
the sensitivity S2 = 2Si. With increasing |fei2| the sen- 
sitivities associated with different channels are canceled 
inside the region [si . . . 82] due to the different signs in 
the charge of the first and the second level and partially 
added outside this region. Given the large level separa- 
tion |fei2| ^ r, the maximum sensitivity drops to the 
sensitivity of the single channel case S2 = Si. 

The effects of the Coulomb interaction in the multi- 
channel case > 1 may be separated into the effect of 
degenerate regions (0 < Vn < N ot Vn = A"), which are 
either partially or totally insensitive to the phase, and the 
effect of the repulsion between normal dot levels. Note 
that the size of the regions with Vn > decreases with in- 
creasing Coulomb interactions between levels, which can 
be seen by comparing Figs. 10(a), 10(c), and 10(e) or 
Figs. 10(g), 10(i), and 10(k). The combination of the 
normal level repulsion and the emergence of the insen- 
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FIG. 10: (Color online) Charge-to-phase sensitivity in the ground state ^2,3,4 in {(p,s) space. All parameters coincide with 
Figs. 8 and 9. In the plots (b) and (h) the sensitivity diverges in {(p,s) = (tt, ±r) and sensitivity maximum (minimum) was 
cut by hands. For all N there are a different colorbars. 



sitive regions leads to "oscillations" in the sensitivity as 
a function of (p and e as presented by the blue and red 
color in Fig. 10. 

The sensitivity as shown in Figs. 10(b), 10(d), 10(f), 
10(h), 10(j), and 10(1) can either increase or decrease due 
to the CAR. Partially, the CAR leads to a divergence in 
sensitivity at the points {(p^ e) = (tt, ±r) in Figs. 10(b) 
and 10(h). 

The sensitivity of the two parallel SWNTs is shown in 
Figs. 10(m)-10(o). 



IX. CONCLUSION 



In this article we have described an Andreev quantum 
dot with several normal levels/conducting channels in the 
infinite superconducting-gap limit. The scaling of the 
charge and the current with the number of channels is 
the central question of the work. We have introduced a 
recursive scheme for an N channel Hamiltonian and have 
analyzed it numerically /analytically. 
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This approach has ahowed us to specify the degeneracy 
of the ground state, depending on the superconducting 
phase difference and the position of the back gate. Due 
to the Coulomb interaction inside each channel, doubly 
degenerate ground states appear. For the case of the 
intermediate Coulomb interaction between the channels, 
regions with a higher degeneracy are generated. Nev- 
ertheless, when increasing the Coulomb interaction be- 
tween the channels, the doubly degenerate regions "re- 
pel" each other, and regions with higher degeneracy dis- 
appear. The size of the doublet regions decrease due to 
the Coulomb interaction between channels. The crossed 
Andreev reflection smears the borders between regions 
with different degeneracy and can completely destroy the 
regions with higher degeneracy. Also, we have studied 
the magnetic properties of the degenerate ground states 
in the presence of a Zeeman splitting. 

Finally, the charge of the intermediate region and the 
current of the multichannel Andreev dot have been com- 
puted. The interplay between the scaling and interac- 
tion effects has been discussed for some realistic situa- 
tions. The charge-to-phase and current-to-gate voltage 
sensitivities increase with the number of channels but do 
not scale linearly as in the case of independent channels. 



While the sensitivity always increases with the number 
of channels, the Coulomb interaction between channels 
leads to a sensitivity reduction. A similar behavior has 
been detected for the critical current as a function of 
the number of channels. The multichannel device could 
therefore be used as the sensitive magnetic flux detector 
or alternatively as the Josephson transistor. Summing 
up, the charge, the current, and the sensitivity scales lin- 
early in the absence of crossed Andreev reflections and 
Coulomb interaction between channels. The Coulomb in- 
teraction suppresses the phase-dependent part, and CAR 
smears the jumps as a function of phase difference and 
gate voltage. 
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